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Abstract 


Theories  and  computational  procedures  are  reviewed  for  processes 
involving  bond  breaking  and  formation  at  a  solid  surface.  These 
processes  include  reactive  scattering,  recombination,  adsorption 
and  desorption.  The  article  ends  with  a  discussion  of  theoretical 
techniques  for  describing  how  some  of  the  above  processes  are  induced 
or  modified  by  laser  radiation. 
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I.  INTRODUCTION 

Among  the  chemical  reactions  of  interest  to  the  chemical  industry 
and  community  in  general,  the  number  occurring  at  a  gas-solid  interface 
is  far  greater  than  that  in  the  pure  gas  phase.  However,  the  most 
thorough  analyses  of  reaction  dynamics  at  the  microscopic  level  have 
been  performed  for  gas-phase  processes,  both  experimentally  and  theo¬ 
retically.  In  regard  to  experiments  on  gas-surface  reactions,  there 
are  special  difficulties  such  as  preparing  a  clean,  specific  crystal  face 
which  are  not  present  in  gas-phase  situations.  The  theoretical  problems 
arise  from  the  many-body  aspect  of  surface  processes,  where  one  must 
deal  with  the  motions  of  many  ,nore  nuclei  than  say,  for  a  three-body 
rearrangement  gas-phase  reaction.  Nevertheless,  strong  progress 
during  the  past  five  to  ten  years  has  been  made  in  the  laboratory  and 
on  the  development  of  appropriate  theories  for  understanding  the  detailed, 
microscopic  dynamics  of  reactions  at  a  gas-surface  interface.  Aside  from 
representing  a  major  area  of  basic  research  in  chemical  kinetics  and 
physical  cherni.--. such  studies  are  leading  to  an  understanding  of 
practical  chemical  processes  such  as  associated  with  heterogeneous 
catalysis. 

In  this  review  article  we  shall  address  some  key  advances  in  the 
theory  of  reactions  at  a  gas-solid  interface.  We  shall  regard  the  term 
"reaction"  in  a  broad  sense  to  apply  to  any  process  involving  bond 
breaking  and/or  bond  formation.  We  shall  direct  our  attention  to  four 
general  types  of  processes:  reactive  scattering,  recombination, 
adsorption  and  desorption.  The  first  of  these  is  discussed  in  Sec.  II. 

Two  different  types  of  reactions  are  considered:  (1)  a  gaseous  molecular 


ion  reacts  with  surface  atoms  via  several  steps,  and  (2)  a  gaseous  atom 
undergoes  a  direct  reactive  encounter  with  an  adspecies.  Both  quantum 
mechanical  and  classical  models  are  utilized,  where  the  latter  involves 
the  intergration  of  classical  trajectories  leading  to  the  solution  of 
the  Generalized  Langevin  Equation.  The  second  process  (Sec.  Ill)  is 
actually  a  combination  of  the  other  three,  since  it  involves  adsorption 
of  a  gaseous  molecule  to  result  in  atomic  adspecies  and  then  recombina¬ 
tion  and  desorption  to  form  to  a  gaseous  molecule  again.  Similar  techniques, 
based  on  potential  energy  hypersurfaces  for  the  collision  dynamics,  are  used 
to  describe  these  first  two  processes. 

The  third  process  is  discussed  in  Sec.  IV,  and  here  the  role  of  excited 
electronic  states  and  the  coupling  of  nuclear  and  electronic  degrees  of 
freedom  are  included  in  a  description  of  the  adsorption  event.  The  fourth 
process  where  an  adsorbed  species  leaves  the  surface  is  considered  in 
terms  of  single-phonon  and  also  multiphonon  transitions  within  the  solid. 

Finally,  recent  theoretical  advances  in  describing  how  laser 
radiation  induces  or  modifies  surface  rate  processes  are  reviewed. 

Models  for  laser-induced  desorption  are  presented,  and  some  ideas  on 
how  the  excitation  of  surface  states  in  semiconductors  and  metals  affects 
charge  exchange  processes  such  as  ion  neutralization  are  discussed. 
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II.  REACTIVE  SCATTERING 

A.  Formation  of  Binary  Compounds  M_.X^:  Bombarding  Clean  Surfaces  of 
Pure  Metal  M  with  Beams  of  Diatomic  Homonuclear  Cations  Xq 

A  simple  quantum  mechanical  model1  has  been  developed  to  study  a 


special  class  of  gas-surface  reactions,  namely 

oM(s)  +  |  X+(g)  -  MaX6(s),  (1) 

where  X2  is  a  diatomic  homonuclear  cation  and  M  a  metal.  The  overall 
reaction  is  modelled  by  a  four-step  process 

X2(g)  +  e'(s)  -  X2(g)  (2) 

X2(g)  -  2X(g)  (3) 

x(g)  -  x(s)  (4) 

oM(s)  +  sx(s)  -  M  X.(s),  (5) 

CX  p 


where  g  and  s  denote  the  gaseous  and  solid  phases,  respectively.  The 
development  of  the  model  is  based  on  the  physical  interpretation  of 
experimental  observations.  The  first  step  describes  the  incoming 
cations  neutralized  via  various  possible  electron  transfer  processes. 

The  second  step  represents  impact  dissociation  of  X2  molecules.  The 
third  step  describes  the  penetration  of  the  nascent  X  atom  into  the 
metal,  its  de-excitation  and  thermal izati on.  Finally,  the  bond  of 
interest  is  formed. 

1.  Neutralization  of  X^ 

The  neutralization  process  is  treated  by  a  simplified  quantum 

2 

mechanical  procedure  involving  resonance  (RN)  and  Auger  (AN)  processes: 
(RN)  x2(u\  v+,  J+)  +  e’(s,  e,  kj)  -  X2(v,  V,  0) 


(6) 
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(AN)  X2(y+,  v+,  J+)  +  e  (s,  e,  kj)  +  e*(s,  c2k2)  -*■  X2(y,  v,  J)  + 

®  ^ ’  E3^3  ^  ( 7 ) 

Here,  e~(s,  c^.)  i  =  1,  2  denotes  an  electron  in  the  metallic  valence 
band  with  energy  e^  and  momentum  k.  and  e’(c,  E^k-j)  denotes  an  electron  in 
the  conduction  band  or  vacuum  continuum  of  states  characterized  by  energy 
e^  and  momentum  k^.  (y+,  v+,  J+)  and  (y,  v,  J)  are  the  electronic,  vibra¬ 

tional  and  rotational  states  of  the  cations  and  molecules,  respectively. 

The  transition  rates  for  these  processes  are  estimated  by  the  Fermi  golden 
3 

rule.  With  the  assumptions  that  only  the  electronic  component  of  the 
cation-metal  long-range  interaction  is  critical  for  these  rates,  that  the 
electronic,  vibrational  and  rotational  motions  of  X^  and  X2  are  separable 
even  at  small  distances  away  from  the  surface,  and  that  the  lattice  vibra¬ 
tions  do  not  participate  in  the  neutralization  processes,  the  transition 
rate  is  given  by 

Ar[  i  -*■  f;  £  (t )  ]  =  -^-|<y+,  Ejkj  |0r(  t)  |  y>j2|  <y+,  v+J+|y,  vJ>|2 

+  +  +  (8) 

6[E(y  v  J  )  -  E(yvJ),  £j] 

Afl[i  **■  f »  i(t)]s  ^  |  <y  ,  e2^2^a^£^  I E3^3>l 

| <y+v+J+jyvJ> |26[E(y+v+J+) -  E(yvJ),  +  c2  -  e^]  (9) 

Here,  6^  _(0  denotes  the  electronic  interaction  operators  responsible 
for  the  resonance  neutralization  and  Auger  neutralization,  respectively. 
i  *  £(t)  is  the  distance  between  the  cation  and  the  surface,  E(y+v+J+) 
and  E(yvJ)  are  the  energies  of  the  cationic  and  molecular  rovibronic 
states,  and  6  is  a  delta  function.  According  to  the  assumption  of 
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separability,  the  square  modulus  of  the  overlap  integral  <y+v+j+|yvJ> 
is  given  by  Franck  Condon  factor  S(y+v+,  yv)  and  Honl-London  factor 
F ( vi  v  J  ,  yvJ):^ 


|<:v+vt)  |uvO>  *  S(y+v+,  yv)  F(y+v+J+,  yvJ), 


(10) 


Integrating  the  transition  rate  with  time  yields  the  transition  probability, 
P^(i  f),  2  s  r,a.  The  probabilities  P2(y+v+J+-*-  yvJ)  are  the  sums  of  the 
transition  probabilities  P  (i  -*•  f),  with  the  summation  being  extended  over 
all  metal  states  whose  energies  are  comparable  with  the  energy  conservation 


ej  *  E(y+vV)  -  E(yvJ)  for  RN  (11) 

n  .  ez  -  e3  -  E(y+v+J+)  -  E(yvJ)  for  AN  (12) 

Combining  Eqs.  (8)  -  (12),  the  probabilities  Pz(y+v+J+  -*•  yvJ)  can  be 
written  as 


P_(y+v+J+,  yvJ)  *  tt  S(y+v+,  yv)  F(y+v+J+,  yvJ) 


e1k1iOr(Oiu>|2 


fi[E(y+v+J+)  -  E(y ,  vJ),  Cj], 


(13) 


P  (y+v+J+,  yvJ )  *  -tt  S(y+v+,  yv)  F(y+v+J+,  yvJ) 


/Gfde,  JCfde-J  d tJ  dkj  dkJ  dk, 
Ub  \itb  3\^  3 


fcb  *  ‘b 

rV 
b0 


/  fdti<y+,  c1k1C2k2|68(i(t)) |y,  c3k3> J2 
to 
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X  6 [ E ( u  ,  V+,  J  )  -  E(yvJ),  +  e2  "  e3-^' 


In  Eqs.  (13)  and  (14),  the  highest  and  lowest  energies  of  the  metal 


valence  band  are  denoted  by  e^  and  cb,  and  /(E.jdk^  denotes  the  integra¬ 


tion  over  all  the  possible  directions  of  the  momentum  k^  associated  with 


the  energy  level  ty  All  the  dynamic  information  is  contained  in  the 
function  £(t).  Both  processes  are  assumed  to  occur  within  the  range 
£(tg)  -  t(t^).  Under  the  experimental  condition  that  the  incoming 
kinetic  energy  is  small,  $  30  eV,  and  that  the  time  of  flight  throe 
the  region  of  long-range  interaction  is  long  enough  whereby  the  con 
neutralization  (RN  +  AN)  is  nearly  100*  efficient,  a  knowledge  of  £  ' 
is  unnecessary,  and  the  integration  over  time  can  be  performed  by  the 
mean-value  theorem,  namely, 


id 


J  fdt  Ar^ectronic(£(t))  =  Are^ectronic(I)  x  |t(J  -  tfl 


where  £(tf)  <  £  £(tQ),  and  Are^ectronic  are  electronic  components  of 


the  transition  rate.  Finally,  A  e]ectronic(T)  are  replaced  by  their 

r ,  a 


"state  averages",  denoted  by  A„  _.  Introduction  of  these  approxima- 

r ,  a 


tions  into  Eqs.  (13)  and  (14)  yields  the  following  expression  for 


P  (y+v+J+  -*■  yvJ): 

i  ,  a 


Pr,.<“Vj*  *  ‘  T  1‘0  -  tfl7Sr..N(“+>Nl“>  *  ar.»(T> 


S(y+v+,  yv )  F(y+v+J+,  yvJ )  , 


Here,  N(y+)  and  N(y)  are  the  multiplicities  of  the  y^th  cationic  and 


+  +.+, 


y-th  moleuclar  states,  and  T  equals  E(yvJ)  -  E(y  v  J  )•  The  function 


(14) 


(15) 


(16) 
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a(T)  is  defined  as  follows: 

'  »o 


er(T)  =  ps(T) 


(17) 


2e*-T 


5a(T)  =  /  f/T,d£l  /  f/T  ,de2  /  f  dc3  ps(el)psU2)pc(£3)’ 


cjd) 


^2  ( ) 


(18) 


where  p$(c)  is  the  density  of  electronic  states  in  the  metallic  valence 
band  and  pcU)  is  the  density  of  states  in  the  conduction  band  and  vacuum 
continuum.  The  integration  limit  is  defined  as  the  following:  the 
minimum  energy  for  an  electron  either  in  the  conduction  band  or  in 
the  vacuum  continuum  is  the  Fermi  energy  e^,  and  the  maximum  energy  is 
given  by  Eq.  (12), 


Emax  _  e^(Se^)  +  e2(=Cf)  -  T.  (19) 

Similarly,  the  limits  for  the  e2  and  integrals  are  given  by 

-  max(e.p  +  T  -  e^»  e^)  <  e2  <  e^  (20) 

and 

e^(T)  *  max(T ,  Eb)  <  £j  <  Ef,  (21) 

respectively.  Information  about  p$(e)  and  pc(e)  can  usually  be  obtained 

e . 

experimentally  from  Auger  electron  spectroscopy*  or  theoretically  such 

5 

as  from  the  free-electron  model.  The  neutralization  probabilities  are 

then  calculable  quantities,  provided  the  "state-averaged"  matrix  elements 

T  ,  are  known. 
r,a 


2.  Phenomenological  Treatment 

The  probability  of  X2  dissociation  is  assumed  to  depend  on  the 
kinetic  energy  Ek  and  the  dissociation  energy  D  =  D(yvJ),  but  insensi¬ 
tive  to  other  specifics  of  the  (yvJ)  state.  A  statistical  model  of 
Kassel^  can  be  applied  to  the  dissociation  probability, 

P(Ekj  D)  =0  for  Ek  <  D 

P(Ek,  D)  =  1  -  (^-)p_1  for  Ek  >  D, 

in  which  p  is  the  number  of  normal  vibrational  modes  of  the  complex. 
This  model  is  justified  if  the  diatom  forms  some  kind  of  complex  with 
the  metal.  Let  R(X2,  Ek)  denote  the  dose  of  X2  cations  in  the  molecul 
beam  at  the  kinetic  energy  Ek  and  W(y+,  v+,  J+)  present  in  the  beam  at 
the  distance  £(tQ)  from  the  surface.  Then  R(X,  Ek),  the  dose  of  atoms 
received  by  the  metal  surface  is  given  by 

R(X,  Ek)  =  2Pt  fdD  P(Ek,  D)  R(X2,  ER,  D), 

R(X2,  Er,  0)  =  [Qa(Ek,  0)  +  Qr(Ek,  D)]  R(xJ,  Efc) 

Qr  , ( E.  ,  D)  =  l  l  W(u+v+J+)  P  >+vV  -*•  yvJ)  6[D(yvJ)-D] 

I  ,G  K  .  .  .  '  *» 

u  v  J  yvJ 

and  Pt  is  the  j>robabil ity  that  an  X  atom  produced  by  Reaction  (3) 
will  reach  the  surface  and  be  de-excited  and  thermal ized.  The  factor 
of  2  in  Eq.  (24)  arises  from  the  stoichiometric  factor  of  Eq.  (3). 
Also,  R(X2,  Ek,  D)  is  assumed  to  be  continuous  with  respect  to  D. 
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For  the  case  of  low  initial  impact  energy,  the  product  is 
formed  only  within  the  first  few  monolayers,  so  that  analytical  measure¬ 
ments  can  be  assumed  to  be  able  to  sample  the  full  thickness  of  the 
layer.  Disregarding  attenuation  of  the  measurement  by  outer  layers, 
the  intensity  I(Ek,  R)  is  then  proportional  to  the  number  of  stoichio¬ 
metric  units  formed, 


I(Ek,  R)  =  Kn(R) . 

Here,  n(R)  is  the  number  of  M  Xr  units  formed  within  the  reaction 

Q  6 

volume  (defined  by  the  surface  layer  thickness  and  area  A  of  the 
beam)  by  a  dose  R(X,  E^)  of  X  atoms  received  by  the  surface  area 
A  from  Process  (4).  The  proportionality  constant  K  contains  the 
relevant  instrumental  factors  and  probability  function  of  the 
analytical  measurement. 

The  relation  of  n(R)  to  R( ,  X)  can  be  obtained  by  solving 
the  differential  equation 

dn(R)  *  [n (M )  -  an(R)]  o(X)/A  dR, 


where  n(M)  is  the  number  of  target  atoms  M  in  the  reaction  volume  of  the 
clean  surface,  and  o(X)  is  the  reaction  cross  section  of  X  atoms  with 
lattice  atoms  M.  Rearranging  Eq.  (28),  we  can  write 


d[n(M)  -  an(R)]  =  [n(M)  -  an(R)]  [-  -^-jdR, 


and  upon  direct  integration  we  obtain 


s  dM 


(1  -  exp(-ao(X)R/A)). 


For  low  doses  of  X^  ions  and  low  impact  energies,  we  have  o(X)  R(X)/A  <<  1 
Under  this  condition  we  can  expand  Eq.  (30)  and  safely  retain  only  the 
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the  first  two  terms  to  obtain 
n(R)  =  n(M)o(X)R(X)/A. 

Substituting  Eqs.  (24)  -  (26)  and  (31)  into  (27),  we  have 
I(Ek,  R)  =  CCaGa(Ek)  +  CrGr(Ek)]  R(X+,  ER) 

where  C =  Kn(M)c(X)Pt^z  z  =  r,  a 

and 

Gz(Ek)  =  /  kP(E  ,  D)  Qz(Ek,  D)dD. 

0  k 


(31) 

(32) 

(33) 

(34) 


3.  Remarks 

A  simple  quantum  mechanical  model  is  constructed  to  provide  quali¬ 
tative  understanding  of  the  detailed  steps  involved  in  the  reaction  of 
low-energy  (<  30  eV)  homonuclear  diatomic  ion  beams  with  metal  surfaces. 
In  this  model,  all  the  dynamical  information  is  contained  in  the  two 
adjustable  parameters  Ca  and  C  [see  Eq.  (33)].  The  theoretical  cal- 
culations  involve  the  evaluation  of  the  band  structure  of  the  metal, 
the  Franck-Condon  factors  and  the  Honl -London  factors  of  cations  and 
molecules.  In  addition,  it  has  been  shown  that  the  electronic  structure 
of  the  metal,  cations  and  molecules  are  initial  factors  in  determining 
low-energy  nitridation  of  a  metal  surface  by 

The  present  model  neglects  the  Auger  de-excitation  process  of 
the  projectiles  which  are  neutralized  by  the  resonant  process.  Never¬ 
theless,  since  both  of  the  electronic  matrix  elements  of  the  Auger  and 
the  resonant  processes  are  left  as  adjustable  parameters,  this  assumption 


does  not  affect  the  model.  On  the  other  hand,  the  neglect  of  the 
electronic  interaction  between  the  metal  surface  and  the  cations/ 
molecules,  which  alters  the  equilibrium  distances  of  the  cations 
and  the  molecules  and  the  electronic  and  vibrational  energy  levels, 
might  introduce  some  error  in  the  computation  of  the  Franck  Condon 
and  Honl -London  factors. 


B.  El ey-Ri deal -Type  Reaction:  0(g)  +  C-Pt(lll) 

Bombarding  oxygen  atoms  on  a  carbon-platinum  adspecies-surface 

8 

system  has  been  studied  by  Tully  with  a  stochastic  trajectroy  technique, 

where  only  direct  (Eley-Rideal -type)  reactive  encounters  are  simulated. 

The  interaction  potential  is  taken  to  be  the  London-Erying-Polyani-Sato 
0 

(LEPS)  form  with  empirical  parameters.  The  potential  incorporates 
experimental  information  about  the  gas-phase  CO  molecule  and  about 

C,  0  and  CO  absorbed  on  Pt(lll).  A  brief  account  of  the  stochastic 
formalism  is  given  below. 

1.  Generalized  Lanqevin  Equation 

In  the  situation  that  localized  reaction  events  take  place  in  the 
midst  of  very  large  assemblies  of  atoms,  the  stochastic  trajectory 
technique  can  be  used  to  eliminate  the  vast  majority  of  "uninteresting" 
atoms  and  to  focus  on  the  local  region  of  action.  The  objective  of  the 
stochastic  trajectory  approach  is  to  accurately  describe  the  flow  of 
energy  into  and  out  of  the  local  region  of  action  without  explicit 
inclusion  of  huge  numbers  of  uninteresting  atoms. 

Q 

The  method  discussed  here,  as  used  recently  by  Tully,  is  based 
on  the  pioneering  work  of  Adelman  and  Doll.*®  There  are  two  basic 
assumptions  in  the  method.  The  first  is  that  the  solid  is  harmonic, 
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such  that  the  mass-weighted  coordinates  of  the  atoms  in  the  solid 
satisfy  the  equation 


"i(t)  ’ 

or  in  matrix  notation, 
u(t)  =  -  fi2u(t). 


where 


(35) 


(36) 


(37) 


m.  is  the  mass  of  the  atom  in  the  solid,  and  K . ,  is  the  force  constant. 

I  I J 

The  second  assumption  is  that  the  forces  experienced  by  the  gas 
atoms  depend  on  the  instantaneous  positions  of  only  a  small  number,  n, 
of  local  surface  atoms.  We  designate  these  as  the  "primary  lattice 
atoms".  The  remaining  N-n  solid  atoms  serve  as  a  heat  bath  for  the 
incident  gas  atom  and  the  primary  lattice  atoms.  These  are  "secondary 
lattice  atoms".  One  can  introduce  the  projection  operators  P  and  Q 
which,  respectively,  project  onto  the  primary  lattice  and  the  secondary 
lattice  atoms.  The  operators  have  the  usual  projection  properties 


P  +  Q  *1;  P2  =  P;  Q2  =  Q;  PQ  =  0.  (38) 

Using  these  properties,  we  can  partition  Eq.  (35)  to  yield 

ijp(t)  *  -Q  pp.p(^)  ”  -pQ^Q^^  (39) 

UqU)  s  -?‘QpUp(t)  -  n2QuQ(t),  (40) 

where 

up(t)  =  Pu(t),  fipQ  -  P(?2Q»  etc.  (41) 
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In  the  presence  of  the  gas-phase  atoms,  the  equations  of  motion 
of  the  gas  atoms,  primary  and  secondary  atoms  are  given  by 


x(t)  =  Fx[x(t),  Up ( t ) ,  uQ] 

Up(t)  ~  ~  GppUp(t)  -  fipgUg(t)  +  fp[x(t),  Up(t),  Uq] 


=  "  ?Qp¥p(t)  ■ 


where  x(t)  are  the  mass-weighted  coordinates  of  the  gas  atoms.  and  Fp 
are  mass-weighted  forces,  derived  from  the  gas-surface  interaction 
potential  U  by 


Fx  [x (t) ,  Up ( t ) ,  Uq]  ■  -  jf:  Cx(t),  up(t),  Uq] 

Fp  [x(t  )  ,  Up(t),  Uq]  =  -  [x(t),  Up(t),  Uq]. 

Note  that  the  forces  only  depend  on  the  equilibrium  positions  Uq  of 
the  secondary  lattice  atoms. 

We  can  write  a  formal  solution  of  Eq.  (44)  as 
n  * 

uQ{ t )  *  ujj(t)  -  /  dt'6(t,  f)?ppup(f), 
where  Ug( t)  is  the  homogenous  solution  given  by 

ug(t)  =  cos{r.QQt)u°(0)  +  npQ  sin(gQQt)  u°(0). 


17 


Choosing  the  Dirichlet  boundary  condition,  the  Green's  function 
can  be  expressed  as 

G(t,  t')  =  Ogg  S1  n[?qQ  ( t-t  *  )]. 


Substituting  Eqs.  (48)  and  (49)  into  (47),  we  obtain 


y q ( ^ )  yq(o)  +  ?qq  yq(^) 


'•^ndt'?QQ  Sin^?QQ^t"t'  ^rQP  yp^'^* 


Upon  integration  by  parts,  Eq.  (50)  becomes 


Uq(t)  -  cos(Cqgt)  Uq(o)  +  Cqq  sin(Qqqt)  Uq ( o ) 


?  QQ  ;QP  “pW  +  ?QQ  C0S(?QQt^  pQp  9p^°) 


+  |^dt  ?qq  c®s[f2qq(t-t 1 )]  fiqp  Up(t'). 

Substituting  Eq.  (31)  into  (43)  one  obtains  the  Generalized  Langevin 
Equation  (GLE): 


up(t)  »  -  feff  Up(t)  -  A(t)  Up(0)  -  J  dt'A(t-f)  Up(t') 


R(t)  +  Fp(x(t) ,  Up(t),  Uq), 
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where 


?eff  =  ^PP  ” 

(53) 

*  r.pg  G  qq  C0S(§gQt)  fiqp 

(54) 

R(t)  =  -Hpg  u[J(t). 

(55) 

Since  the  time  dependent  of  Ug(t)  is  influenced  by  the  motion  of  the 
primary  atoms,  the  term  R(t)  behaves  like  an  external  force  in  Eq.  (52). 
Furthermore,  the  initial  conditions  Uq(0)  and  Ug(0)  are  not  known  precisely 
and  are  randomly  distributed  for  a  given  temperature.  Therefore,  R(t) 
can  be  taken  as  a  random  force,  and  because  of  the  random  impulses 
injected  by  R(t),  the  trajectories  are  termed  "stochastic". 

The  above  GLE  is  exact,  but  the  friction  kernel  A(t)  and  the 
random  force  R(t)  are  very  complicated.  Exact  evaluation  of  Eq.  (52) 
would  involve  the  same  labor  as  solving  the  original  enormous  set  of 
Eqs.  (43)  and  (44).  The  objective  of  the  GLE  approach  is  to  approxi¬ 
mate  A(t)  and  R(t)  in  such  a  way  that  they  are  easy  to  compute  and  yet 
adequately  describe  the  effects  of  the  heat  bath  on  the  primary  lattice 
atoms.  First,  R(t)  can  be  represented  by  a  Gaussian  random  force  due 
to  the  assumption  of  the  harmonic  motion  of  the  secondary  lattice  atoms. 
Second,  the  autocorrelation  function  of  R(t)  is  related  to  the  friction 

kernel  A(t)  by  means  of  the  "second  fluctuation-dissipation  theorem,"11 
* 

<R( t )  R+(0)>  *  kBTA(t),  (56) 

where  T  is  the  temperature  of  the  lattice  and  kg  is  the  Boltzmann's 
constant^  Qualitatively,  it  states  that  the  energy  dissipated  from 
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the  primary  zone  by  the  friction  must  balance  the  energy  introduced 
by  the  fluctuating  force  in  order  to  maintain  a  certain  temperature 
T.  Eq.  (56)  provides  us  with  a  prescription  for  constructing  R(t) 
as  long  as  A(t)  has  been  determined;  thus  simulations  can  be  per¬ 
formed  at  any  desired  temperature,  no  matter  how  small  the  primary 
zone.  In  the  actual  computation,  the  friction  kernel  A ( t )  is  usually 
modelled  phenomenologically;  more  precisely,  a  simple  expression  is 
employed  for  A ( t )  with  some  adjustable  parameters  such  that  the  surface 
and  the  bulk  properties  can  be  best  reproduced,  for  example,  £(t)  can 
be  expressed  as  the  solution  of  the  damped  harmonic-oscillator  differential 

equation,  given  by 

A(t)  =  A  exp(->t  +  i  /? -  ?t, , 

where  y..  is  the  damping  coefficient  and  is  the  harmonic  frequency. 
Thus,  the  elements  of  the  three  matrices  A,  y  and  u  are  chosen  to  best 
reproduce  known  experimental  or  theoretical  information  about  the  surface 
vibrations  of  a  lattice. 

2.  Potential  Energy  Hypersurface 

The  interaction  potential  is  the  most  uncertain  and  also  the  most 

crucial  Ingredient  to  this  study.  An  empirical  interaction  potential, 

12 

as  mentioned  above,  of  the  LEPS  form,  based  on  experimental  informa¬ 
tion  about  the  isolated  CO  molecule  and  about  C,  0  and  CO  adsorbed  on 
Pt(lll),  can  be  employed.  However,  due  to  limited  experimental  infor¬ 
mation,  there  is  a  large  uncertainty  associated  with  this  procedure. 
Nevertheless,  the  interaction  potential  can  be  varied  by  adjustable 


(57) 
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parameters  such  that  gross  qualitative  features  of  the  reaction  dynamics 
are  elucidated. 

The  main  interest  in  the  present  reaction  is  the  energy  disposal 
of  the  exothermicity.  The  binding  energy  of  carbon  atom  to  Pt(lll), 
which  directly  affects  the  exothermicity  of  the  reaction,  is  probably 
the  most  uncertain  feature  of  the  interaction  potential.  Reasonable 
values  for  the  C-Pt(lll)  binding  energy  range  from  4.0  to  7.5  eV, 
corresponds  to  reaction  exothermicities  of  7.6  to  4.1  eV.  For  this 
reason,  the  empirical  parameters  of  the  interaction  potential  are  varied 
such  that  the  exothemicity  changes  substantially  from  6.1  to  3.8  eV, 
where  it  is  seen  that  gross  features  of  the  reaction  dynamics  have 
little  dependence  on  the  magnitude  of  the  exothermicity. 

3.  Results  and  Remarks 

The  reaction  of  the  gas-phase  oxygen  atom  with  carbon  absorbed  on 

platinum(lll)  has  been  simulated  in  the  following  way:  oxygen  is 

directed  normal  to  the  surface  with  the  two  fixed  initial  kinetic 

© 

energies,  2  and  10  kcal/mole,  at  an  initial  distance  of  12  A  away 
from  the  surface.  Two  different  surface  temperatures,  0  and  500°  K, 
are  examined.  In  all  cases  studied,  the  reaction  probability  of  forming 
a  CO  molecule  is  very  high.  This  is  no  surprise,  since  there  is  no 
potential  barrier  in  the  interaction  potential.  Furthermore,  most  CO 
molecules  escape  from  the  surface  within  a  few  picoseconds.  A  large 
percentage  of  the  liberated  energy  is  disposited  in  the  gas-phase  CO 
molecule,  partitioned  among  translations,  vibrational,  and  rotational 
degrees  of  freedom.  The  rotational  motion  is  formed  to  be  anisotropic, 
i.e.,  molecules  tend  to  tumble  rather  than  rotate  in  the  plane  of  the 
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surface.  The  angular  distribution  is  broad  and  is  biased  toward  the 

normal  direction.  Significant  alterations  of  the  interaction  potential 

can  result  in  a  change  of  the  partition  of  the  energy  among  translation, 

vibration  and  rotation  of  the  CO  molecules. 

In  the  present  study,  the  energy  has  been  assumed  to  flow  to  and 

from  vibrational  modes  of  the  solid.  However,  contributions  from 

electronic  excitation  and  de-excitation  within  the  solid  have  not 

been  accounted  for.  This  may  introduce  substantial  error  in  the 

present  reaction,  since  platinum  is  a  metal  possessing  conduction 

electrons  that  can  be  easily  excited.  Also,  possible  quantum  effects 

13 

such  as  vibrational  adiabaticity,  which  is  important  in  gas-phase 
reactions,  have  not  been  considered. 


III.  RECOMBINATION 


Atomic  recombination  dynamics  of  the  ^  +  W(OOl)  system  have  been 
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investigated  by  McCreery  and  Wolken  with  a  model  potential  of  the 
LEPS  form.  The  process  involves:  (i)  initial  adsorption  of  the  hydrogen 
as  atoms  on  the  surface,  (ii)  possible  equilibration,  and  (iii)  recom¬ 
bination  and  description  of  the  hydrogen  molecule.  This  sequence  is 
known  as  the  Langmuir-Hinshelwood  process.  The  model  assumes  that  the 
surface  is  rigid,  which  means  that  there  is  no  energy  transfer  between 
the  atom  and  lattice.  The  assumption  should  be  valid  when  the  absorption 
time  and  the  recombination  and  description  times  are  much  shorter  than 
the  absorbate-surface  vibrational  relaxation  time. 

In  this  section,  we  will  consider  the  following  issues:  (i)  the 
formation  of  the  rigid  surface  potential,  (ii)  the  extention  to  the  moving 
surface  model  with  the  lattice  dynamics  described  by  the  Einstein  model; 
and  (iii)  comparison  of  the  Einstein  model  with  the  GLE  for  different 
values  of  macroscopic  parameters. 

A.  Rigid  Surface  Model  Potential 

The  LEPS  potential  is  derived  from  a  valence  bond  treatment  for  a 
system  of  three  one-electron  atoms.  For  the  problem  of  a  diatom  AB 
colliding  with  a  solid  surface,  there  are  also  three  two-body  interactions 
of  interest:  atom  A-surface,  atom  B-surface  and  atom  A-atom  B.  Due  to 
the  infinite  extent  of  the  solid,  the  interbody  separations  are  subjected 
to  fewer  geometrical  constraints  than  the  three-atom  case,  and  a  new 
asymptotic  limit  is  introduced  where  both  atoms  are  near  the  surface  but 
the  A-B  interatomic  distance  is  large. 
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The  energy  expression  for  the  four-electron  four-orbital  valence 
bond  structure  is  given  by  ^ 

E  =  Qu+Qj  +  Q  +Qi.J  +  Q^+QJ 
vab  ycd  yac  vbd  vcb  wad 


-  {*r((ci  "  8)^  +  (y  -  &)^  +  (s  -  y ) ^ ] ) 2 


(58) 


where  is  the  Coulombic  energy  for  the  two-electron  system  i-j,  and 


a  =  a-|  +  <*2 

=  (ab). 

a2  =  (cd) 

e  “  b1  +  e2 

£  1  =  (ad). 

B2  =  (be) 

(59) 

'  =  Y1  +  y2 

^  =  (ac), 

y2  =  (bd) 

are  exchange  integrals  for  the  two  centers  indicated.  More  precisely, 

fab)  =  <abcd|H'bacd>,  etc.  (60) 

For  the  case  of  a  diatomic  molecule  interacting  with  a  solid  surface,  the 
valence  bond  energy  formula  is  modified  such  that  the  solid  can  provide 
the  electrons  for  bonding  to  the  incident  atoms  at  any  point  (or  points) 
on  the  surface.  Therefore,  a  simple  model  of  chemisorption  is  obtained 
by  only  considering  those  structures  describing  a-b,  a-d,  and  b-c  inter¬ 
actions.  This  leads  to  the  energy  expression 

E  ■«,„  +  «,d  *  «bc  -  4I(oi  - 6)2  +  •*»'■  <6’> 

We  follow  the  usual  LEPS  treatment  by  assuming  that  the  atom-atom  and 
atom-surface  interactions  are  adequately  described  by  Morse  potentials. 

Additional  flexibility  is  gained  by  introducing  a  Sato  parameter  a,  which 
is  inserted  in  $ych  a  way  as  to  have  no  effect  when  one  or  more  of  the 
distances  becoipes  large.  Thus  we  obtain  a  diatom-solid  potential  of 
the  form 
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V  =  T?A  {V1  +  V2  +  V3  "  [A1  +  (A2  +  A3)2  '  A1(A2  +  A3)]4}  (62) 

Pi 

where  Vi  =  -5-  [(3+A)  exp(-2&i  ir.-r.Q))  -  (2+6A)  exp(-a.(ri-r.0))]  (63) 

Di 

Ai  =  -4  [(H3A)  exp(-2ai(ri-rio))  +  (6+2A)  exp(-a.  (^-r^))]  (64) 

and  D.j ,  a^  and  r.  are  the  dissociation  energy,  Morse  parameter  and 
equilibrium  distance  for  the  i-th  two-body  interaction.  The  power  of  the 
LEPS  method  lies  not  so  much  in  its  quantitative  description  of  the  po¬ 
tential  as  in  a  reasonable  description  of  the  potential  surface  with 
sufficient  flexibility  to  permit  the  study  of  a  wide  variety  of  shapes 
of  potentials. 

To  describe  the  static  periodicity  of  the  solid  surface,  the 
dissociation  energy  and  the  equilibrium  distance  of  the  H-W  (001)  are 
written  as  functions  of  x  and  y,  which  are  Cartesion  coordinates  parallel 
to  the  plane  of  the  W  (001)  surface,  given  by 


D(x,y)  =  D0[l  +  6  Q(x,y)] , 

(65) 

rQ(x,y)  =  Zm[l  +  e  P(x,y)j. 

(66) 

Q(x,y)  and  P(x,y)  are  chosen  to  have  the  correct  synwetry  for  the  W  (001) 
surface  and  enough  adjustable  parameters  to  fit  the  Extended  Hiickel 
Molecular  Orbital  (EHM0)  computed  points^  for  the  three  possible  binding 
sites,  1CN,  2CN  and  5CN,  on  the  surface.  The  numbers  1,  2  and  5  correspond 
to  the  number  of  tungsten  surface  atoms  coordinated  to  the  adsorbed  hydrogen 
atom.  Thus,  Q(x,y)  takes  the  simple  form 
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Q(x,y)  =  cos (~— ■)  +  cos(^)  -  A[cos(^p-)  -  lJ[cos(^)  -  1]  ,  (67) 

where  "a"  is  the  1CN  separation.  Given  that  DQ  is  the  dissociation  energy 
at  the  2CN  site,  6  is  then  determined  by  the  binding  energy  at  a  1CN  site 
and  Dq,  and  A  is  determined  by  the  binding  energy  at  a  5CN  site  and  6. 

P(x,y)  takes  a  form  similar  to  Eq.  (67)  but  with  constants  determined 
from  the  equilibrium  heights  above  the  1CN,  2CN  and  5CN  positions. 

B.  Moving  Surface  Model  Potential 

The  rigid  surface  potential,  V.^,  described  above  can  be  generalized 
to  a  moving  surface  model  by  introducing  two  correction  terms: 

(1)  VR,  the  restoring  force  on  each  atom  of  the  solid,  tending  to 
return  it  to  its  lattice  site; 

(2)  VC0RR,  to  account  for  the  change  in  the  gas-surface  interaction 
when  the  solid  atom  is  displaced  from  its  lattice  site. 

Thus,  we  have 

VTotal  =  VLEPS  +  VR  +  VC0RR'  (68' 


There  are  two  basic  assumptions  in  this  model.  First,  there  is  only  a  small 
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number  of  moving  surface  atoms.  Second,  the  Einstein  model  is  used  to 
describe  the  restoring  force  of  the  moving  surface  atoms,  whereby  VR  is 
represented  by  a  harmonic  potential  binding  an  atom  of  the  solid  to  its 
(fixed)  lattice  site.  On  the  other  hand,  Morse  potentials  for  VcqRR 
connecting  each  gas  atom  with  each  surface  atom  are  used  to  account  for 
the  change  in  the  gas-solid  potential  due  to  displacement  of  the  surface 
atom  from  its  lattice  site.  However,  V^p^  already  contains  the  inter¬ 
action  of  each  gas  atom  with  each  lattice  site.  To  avoid  including  it  twice, 
this  lattice  site  interaction  must  be  substracted  from  Vj.qRR,  so  that 
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VCORR"  ^  l  CW(Rg-s)  *  W<Rg-l)]*  (69) 

g  S  3  3 

where  N  is  the  number  of  gas  atoms,  N  is  the  number  of  moving  surface 
y  5 

atoms,  and  R  is  the  distance  between  the  gas  atom  and  the  lattice  site 
of  the  solid  atom.  As  required,  if  the  moving  surface  atom  occupies  its 
lattice  site,  i.e.,  R$_g  =  Rs_-j»  then  V^qRR  =  0.  The  total  gas-solid 
potential  for  the  moving  surface  case  is  then 


where  VR  is  the  harmonic  restoring  force  for  the  moving  surface  atoms. 
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McCreery  and  Wolken  have  used  this  moving  surface  model  to  examine 
the  validity  of  the  rigid  surface  approximation  for  the  recombination 
dynamics  of  H2  +  W  (001).  Out  of  100  trajectories,  only  a  small  percentage 
of  the  total  energy  available  to  the  product  molecules  is  transferred 
between  the  gas  atoms  and  the  surface  atoms.  Furthermore,  the  energy 
distribution  of  the  escaping  molecules  is  very  similar  for  both  the  rigid 
surface  and  the  moving  surface.  However,  the  amount  of  recombination  for 
the  moving  surface  is  about  twice  as  large  as  for  the  rigid  surface.  This 
suggests  the  inadequacy  of  the  rigid  surface  representation  for  describing 
the  recombination  dynamics.  An  unanswered  question  is  the  accuracy  of  the 
Einstein  model.  To  fully  address  the  problem,  a  complete  lattice  dynamics 
calculation,  such  as  with  the  GLE,  is  required  to  justify  the  Einstein 
model  and  the  rigid  surface  assumption. 
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C.  Comparison  of  the  Einstein  Model  with  the  Generalized  Lanqevin 
Equation  for  Different  Macroscopic  Parameters 

29 

Diebold  and  Wolken  studied  the  energy  accomodation  and  dissociative 

adsorption  of  diatomic-surface  scattering  on  a  LEPS-type  potential  with  two 

different  methods.  One  is  the  GLE  approach  described  in  Sec.  II,  but  with 

an  additional  approximation:  the  friction  kernel  A(t)  is  approximated  by  a 

single-term  damped  cosine  function  obtained  by  a  numerical  fit  to  a  modified 

bulk  Debye  model  of  the  solid.  The  other  method  is  the  Einstein  model 

approach,  which  simply  replaces  the  friction  kernel  and  the  random  force 

by  a  single  oscillator.  The  interaction  potential  was  taken  to  be  the  form 

21 

of  Gelbs  and  Cardillo  describing  the  interaction  of  H2  with  Cu  (001). 

The  hydrogen-copper  potential  parameters  were  chosen  to  be  those  of 
22 

Gregory  et  al . 

Two  different  collison  configurations  -  one  that  leads  to  dissociative 
adsorption  and  one  that  does  not  -  were  examined  for  different  solid  tem¬ 
peratures,  Debye  temperatures  and  atomic  masses.  For  the  latter  collision 
configuration,  one  restricts  the  bond  axis  of  the  molecule  to  be  perpendicular 
to  the  plane  of  the  surface  and  positions  the  incoming  molecule  directly  over 
the  struck  solid  atom.  Here  the  energy  exchange  is  directly  analyzed  without 
any  complication  due  to  multiple  collisions  with  the  surface  and  diatom 
rotational  dynamics.  The  internal  energy  of  the  diatom  was  chosen  to  be 
the  same  as  that  of  the  ground  vibrational  and  rotational  state  of  the 
hydrogen  molecule.  For  a  fictitions  hydrogenic  molecule  composed  of  atoms 
twenty  times  as  massive  as  a  hydrogen  atom,  the  amount  of  energy  transfer 
is  seen  to  decrease  slowly  with  increasing  solid  temperature  for  both  the 
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single-oscillator  and  GLE  approaches.  However,  the  results  of  the  GLE 
approach  are  generally  found  to  be  twice  as  large  as  those  of  the  single- 
oscillator  approach.  Increasing  the  Debye  temperature,  at  a  fixed  solid 
temperature,  the  energy  exchange  is  seen  to  decrease.  As  before,  the 
results  of  the  single  oscillator  are  too  small  compared  to  those  of  the 
GLE.  For  a  300°K  solid  with  the  copper  Debye  temperature,  energy  transfer 
of  the  hydrogen  molecule  and  fictitious  molecules  with  atoms  ten  times 
and  twenty  times  the  mass  of  hydrogen  were  examined.  The  results  of  the 
hydrogen  molecule  are  very  similar  for  both  the  GLE  and  single-oscillator 
approachs,  but  the  agreement  starts  to  deviate  as  the  atomic  mass  increases. 

For  the  former  collision  configuration  leading  to  dissociative 
adsorption,  the  sticking  probabilities  were  computed  as  functions  of  those 
macroscopic  parameters.  In  all  cases  studied,  the  single-oscillator  model 
shows  the  same  trends  as  the  GLE  results,  but  with  much  smaller  sticking 
probabilities.  There  are  two  exceptions:  the  single-oscillator  model 
converges  to  the  GLE  results  as  the  atomic  mass  decreases  or  the  Debye 
temperature  increases.  Thus,  one  can  conclude  that  the  single-oscillator 
model  is  a  good  approximation  when  the  atomic  mass  is  small  or  the  incident 
velocity  is  large. 
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IV.  ADSORPTION 

Microscopic,  dynamical  treatments  of  the  adsorptive  process  are 
still  at  an  early  stage  of  development,  and  such  treatments  that  do 
exist  are  restricted  either  to  classical  trajectory  studies  using 
effective  potentials,  or  to  highly  simplified  quantum  or  semi -classical 
models.  What  all  these  treatments  have  in  coninon  is  some  mechanism 
for  dissipation  of  the  excess  kinetic  energy  of  the  incoming  particle, 
and  they  may  be  classified  accordingly.  We  list  here  four  important 
classes  of  dissipation: 

(1)  phonon  excitation 

(2)  polarization  effects 

(3)  dissociation 

(4)  rotational  excitation. 

[Of  these,  the  last  is  not  strictly  dissipative,  in  the  sense  of  an 
irreversible  decay,  but  is  included  for  completeness.]  There  is  very 
little  literature  on  dissociation  as  a  mechanism  for  dissipation,  such 
as  there  is  being  devoted  mainly  to  static  problems.  We  therefore  make 
a  survey  of  work  falling  into  classes  (1),  (2)  and  (4). 

A.  Dissipation  by  Phonon  Excitation 

The  methodology  for  the  most  detailed  calculations  to  date  has  been 

9  10 

the  Generalized  Langevin  Equation,  ’  which  is  reviewed  in  Sec.  II  of 
this  article.  In  other  words,  the  most  quantitative  calculations  have 
been  classical  treatments  using  effective  (dissipative)  potentials.  In 
principle,  nothing  further  needs  to  be  added,  since  in  this  instance  the 
unifying  concept  for  all  the  various  kinds  of  surface  process  is  the 
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classical  trajectory,  the  only  difference  between  adsorption  and 
desorption,  say,  being  that  one  trajectory  ends  where  the  other  begins. 
It  is  only  in  quantum  mechanics  that  different  kinds  of  model  are  used 
to  treat  different  processes. 

A  phenomenological  semi -classical  treatment  has  recently  been 
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given  in  which  the  technology  of  gas-phase  collision  theory  has  been 
transferred  to  treat  adsorption.  The  incoming  particle  is  assumed  to 
move  on  a  repulsive  potential  with  respect  to  the  surface  and  interacts 
via  an  avoided  crossing  with  a  bound  electronic  state.  The  system  can 
therefore  form  resonances,  but  the  adatom  will  eventually  dissociate. 

If  however  the  trajectory  includes  dissipation,  then  after  oscillating 
a  few  times  in  the  bound  state  the  particle  drops  below  the  crossing 
point  and  is  effectively  adsorbed.  The  model  is  a  multiple-pass  Landau 
Zener  problem,  with  dissipation.  In  the  case  of  no  dissioation,  the 
probability  of  returning  out  is  given  by 
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where  is  the  landau-Zener  probability  of  switching  curves  at  each 
pass  of  the  crossing-point.  Clearly,  what  goes  in  must  come  out.  However, 
a  crude  approximation  for  the  dissipation  is  given  by  restricting  the 
number  of  passes  to  a  finite  nuntoer  N.  This  is  calculated  by  assuming 

r 

a  constant  decay  rate  ^  ,  that  is,  a  constant  rate  of  decaying  below  the 

* 

crossing-point.  Then  the  time  taken  to  decay  is  just  t  *  p  ,  and  the 
number  of  passes  that  occur  before  this  happens  is  just 
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N  -  vt  «  ^  ,  (72) 

where  v  is  the  vibrational  frequency  of  the  bound  motion.  It  is  clear 
that  this  phenomenological  treatment  is  quite  limited  unless  values  for 
Pj2  and  r  are  known.  This  would  in  general  be  the  most  difficult  part 
of  the  problem.  However,  there  is  a  more  severe  limitation  on  this 
kind  of  treatment  for  phonon  dissipation,  namely  an  implicit  assumption 
of  weak  coupling.  When  the  coupling  is  strong  it  is  impossible  to 
separate  the  motion  of  the  adatom  from  the  motion  of  the  surface.  Conse¬ 
quently,  a  proper  representation  must  in  general  go  beyond  such  a  simole 
two-body  Landau-Zener  treatment,  in  which  the  many-body  aspects  have 
been  bundled  into  a  phenomenological  width  term.  We  note,  however,  that 
the  above  treatment  is  so  general  that  it  could  be  applied  to  problems 
with  completely  different  mechanisms  of  dissipation. 

B.  Dissipation  by  Excitation  of  Electron-Hole  Pairs 

The  phenomenon  of  dissipation  by  polarization  effects  is  particularly 
appropriate  for  deposition  on  metal  surfaces,  where  such  effects  are  large. 

This  class  would  include  for  example  dissipation  due  to  plasmon  excitation. 
However,  in  the  following  we  restrict  ourselves  to  a  discussion  of  electron- 
hole  pair  formation. 

We  now  consider  the  process  by  which  an  adatom  couples  with  the  bulk 

electrons  of  a  metal  causing  excitation  of  the  substrate.  The  literature 
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on  this  effect  is  very  extensive  and  is  based  mostly  on  the  semi- 

25 

classical  Anderson  model,  although  there  have  recently  been  attempts 

OA  O  Q 

to  introduce  fully  quantum  descriptions.  ’  We  shall  describe  the 
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semi -classical  treatments.  The  overall  process  is  one  of  considerable 
complexity  and  consequently  all  treatments  to  date  have  made  extensive 
use  of  models  and  approximations.  The  essential  idea  is  that  an  atom 
approaching  a  surface  causes  a  time- varying  perturbation  which  can 
promote  electrons  near  the  Fermi  surface  to  excited  states,  thereby 
creating  electron-hole  pairs.  As  a  result,  the  atom  may  lose  an  amount 
of  kinetic  energy  sufficient  to  allow  it  to  stick.  It  is  assumed  through¬ 
out  that  the  atom  maintains  its  character  even  in  the  vicinity  of  the 
surface  and  does  not  form  a  chemical  bond.  The  descriptions  therefore 
are  limited  to  theories  of  physisorDtion.  Moreover,  only  the  attractive 
part  of  the  interaction  potential  contributes  to  electron-hole  pair 

formation,  it  is  being  assumed  that  in  the  repulsive  region  phonon  excita- 
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tion  is  the  dominant  dissipative  mechanism.  It  is  implicitly  assumed 
also  that  the  electrons  interact  with  the  adsorbate  by  means  of  a  local 
potential.  Among  other  things,  this  means  physically  that  individual 
excitation  events  are  independent.  Therefore,  the  total  probability  for 
adsorption  can  be  obtained  by  summing  the  probabilities  of  all  events 
that  can  lead  to  adsorption. 

The  model  is  not  fully  consistent  since  the  trajectory  of  the 
adatom  is  calculated  on  the  assumption  that  no  energy  is  lost.  In  practice, 
however,  this  approximation  is  not  expected  to  be  severe  since  the  strongest 
interaction  occurs  where  the  local  velocity  is  very  much  larger  than  the 
initial  velocity.  Consequently,  the  change  in  velocity  due  to  the  substrate 
excitations  is  expected  to  be  relatively  small.  Finally,  the  semi -classical 
approximation  will  break  down  where  the  classical  trajectory  is  an 
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inappropriate  description  of  the  nuclear  motion,  that  is,  for  light 
adatoms  at  low  collision  energies. 

The  central  quantity  to  be  calculated  is  the  probability  Pe>(c) 
that  an  incoming  particle  with  energy  c'  loses  energy  e  to  the  sub¬ 
strate  during  one  round  trip.  We  assume  that  the  substrate  is  initially 
in  a  state  denoted  by  | > ,  where  e.|  is  the  energy  and  a-j  is  a  set  of 
quantum  nunbers  (such  as  polarization,  spin,  etc.)  which  complete  the 
description.  The  entire  treatment  may  be  given  in  terms  of  standing 
wave  states  since  the  Born  approximation  to  the  electron  scattering  is 

applied.  Furthermore,  let  us  assume  that  the  time  dependence  of  these 
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states  is  given  by  a  coefficient  c  (t).  Then  we  can  write 

e1al 
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where  V12(t)  =  <e2°L2lV(t)  !elal>  ’ 


and  V  is  the  interaction  due  to  the  adatom. 

Vj2(-  00 )  =  0,  we  obtain 
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The  sticking  probability  is  given  by 


Integrating  by  parts,  using 
i(e2-e1)t/tf 


(74) 


00 

s  =  f '  dePe^(e) . 


(75) 


This  means  that  a  particle  with  energy  c'  must  lose  an  amount  of  energy 
equal  to  at  least  e"  in  order  to  stick.  From  here  on  we  suppress  the 
subscript  e',  it  is  being  assumed  that  we  are  considering  only  these 
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particles.  Then  we  can  write 

P(e)  =  l  If  dc1d£2|c£.  (“)  |2f(c1)(l-f(c2))S(e-(G2-c1)),  (76) 

ala2  22 

where  f  is  the  Fermi-Dirac  distribution.  In  other  words,  the  excitation 
can  only  take  place  if  the  initial  state  is  occupied  and  if  the  final 
state  had  been  vacant.  A  very  useful  form  for  the  matrix  element  in 
Eq.  (74)  appears  in  the  Born  approximation.  That  is,  if  the  states  of 
the  electronic  substrate  are  represented  as  ordinary  (time-independent) 
plane  wave  states,  then  we  can  write  it  in  terms  of  the  Born  approximation 
T-matrix 

<e2a2 1 ^( t )  lciai>  =  ft  TU2a2;elal)-  (7~) 

Furthermore,  let  us  assume  that  significant  transitions  occur  only  at  the 
Fermi  energy,  and  that  the  difference  between  and  neqliqible. 

Then 

<C2»2l  V(t)  | e i &i >  'v  ^  T( Ep&2 > ^-p^i ) »  (78) 

where  Cp  is  the  Fermi  energy. 

It  is  always  possible  to  choose  a  representation  which  is  diagonal 
In  a-[  and  o^.  The  matrix  element  in  this  adiabatic  representation  is 
given  by  |^- T(cp0t;cpa).  However,  we  can  write 
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where  6  is  the  elastic  scattering  phase  shift  and  the  approximation 

Ep,a 

is  for  weak  coupling.  Hence 
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(81) 


where 

t  =  e2  -  er 


The  interaction  between  the  adatom  and  substrate  will  be  greatest 

as  the  energy  of  the  atomic  level  approaches  c _.  Then  the  elastic 

scattering  of  electrons  is  dominated  by  the  resonant  interaction  with 

the  one-particle  adatom  state.  Under  these  circumstances  the  phase 
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shift  is  given  by 


6  (t)  =  tan 
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where  F(t)  is  the  width  of  the  adatom  level  and  e  the  real  part  of  its 

O 

energy.  Notice  that  as  the  energy  of  the  level  crosses  the  Fermi  energy, 
the  phase  shift  varies  from  to  -  £  . 

We  made  two  further  assumptions.  First,  r  is  roughly  constant  during 
the  traversal  of  the  Fermi  surface.  We  therefore  set  r( t)  =  r(tQ),  where 
t  is  the  time  of  the  crossing.  Second,  A  is  a  linear  function  of  time 
(this  is  analogous  to  one  of  the  Landau-Zener  approximations).  We  then 
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have 
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where  x  is  the  time  taken  for  the  system  to  pass  the  crossing  point 
(that  is,  for  the  phase  shift  to  change  by  tt).  Notice  that  this  is  a 
Lorentzian  function  of  the  time,  which  is  not  surprising  since  it  has 
been  derived  for  the  case  of  resonance  elastic  scattering.  This  form 
for  6  can  now  be  substituted  back  into  Eq.  (81)  and  the  sticking  proba¬ 
bility  calculated  using  Eq.  (76). 
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A  related  approach  ’  treats  the  electronic  excitations  as  bosons 

and  finds  solutions  to  the  Heisenberg  equations  of  motion  of  the  occupation 

numbers.  This  leads  to  a  simplified  form  for  P(e)  but  appears  at  first 

glance  to  be  an  unwarranted  assumption,  except  in  the  case  of  weak  coupling. 

However,  it  has  been  demonstrated  that  the  form  holds  even  for  strong 
27 

coupling,  provided  the  electron-hole  pair  is  treated  as  a  Tomonaga 
boson. 37 

Realistic  calculations  of  P(e)  show  that  systems  which  undergo  a 

rapid  crossing  of  the  Fermi  level  lend  to  stick  better  than  systems  which 
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make  the  traversal  more  slowly.  At  the  opposite  extreme  are  systems 
which  do  not  undergo  any  crossing,  that  is,  systems  for  which  the  unfilled 
adatom  levels  are  far  removed  from  the  Fermi  surface  at  all  points  on  the 
trajectory.  Thus  the  rare  gas  atoms,  which  have  filled  shells,  stick 
much  less  readily  by  this  mechanism,  and  it  is  likely  that  adsorption  is 
due  to  phonon  dissipation.  In  general,  however,  it  is  necessary  to  consider 
both  electronic  and  phonon  dissipative  mechanisms. 


37 


C.  Rotational  Trapping 

We  consider  a  diatomic  molecule  approaching  a  surface  in  a  potential 
which  can  support  bound  states.  With  no  mechanism  for  dissipation,  it 
will  scatter  out  elastically  after  colliding  with  the  repulsive  part  of 
the  surface-molecule  potential.  However,  during  this  collision  it  may 
undergo  rotational  excitation  in  the  same  way  as  a  football  which  is 
bounced  off  a  wall.  Thus  translational  energy  is  converted  into  rota¬ 
tional  energy,  and  if  the  loss  is  sufficient,  the  molecule  drops  into  a 
bound  state  of  motion. 

38 

R-matrix  calculations  have  been  performed,  for  HD  on  Pt(lll)  and 
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HD  on  Ag(m),  to  assess  the  degree  of  rotational ly  inelastic  scattering. 

The  potentials  used  are  for  rigid  corrugated  surfaces  and  are  of  the  form. 

V(z,e)  =  V0( Z)  [1  +  8P2(cose)J,  (86) 

where  Z  is  the  normal  distance  from  the  surface  to  the  geometric  center 
of  HD,  and  e  is  the  polar  angle  of  the  molecular  axis  from  the  normal  to 
the  surface.  VQ( Z)  is  given  as  a  Morse  potential.  The  results  are  for 
rotationally  elastic  scattering  as  a  function  of  initial  translational 
energy.  The  trapping,  however,  is  dynamic  in  the  sense  that  the  final 
state  is  a  scattering  state. 
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V.  DESORPTION 

When  an  adsorbed  substance  leaves  the  surface  in  a  manner  independent 
of  that  of  its  arrival,  we  may  call  it  desorption  as  opposed  to  scattering. 
As  the  desorption  rate  depends  on  the  surface  coverage  and  activation  energy 
as  well  as  the  temperature,  desorption  may  be  a  useful  means  of  identifying 
and  adjusting  the  surface  population  and  properties  of  adspecies  of  interest 
as  well  as  determining  activation  energy  and  surface  coverages,  th. sorption 
processes  thus  have  fundamental  implications  to  chemical  kinetics  and  mole¬ 
cular  dynamics  on  surfaces. 

For  desorption  to  occur,  energy  needs  to  be  supplied  to  surmount  a 
barrier,  either  due  to  a  van  der  Waals  interaction  for  the  case  of  physi- 
sorption  or  an  actual  bond  for  the  case  of  chemisorption.  The  required 
energy  can  be  provided  by  the  impact  of  ions  or  electrons,  by  electro¬ 
magnetic  radiation  such  as  a  laser  (photodesorption),  or  by  thermal 
heating  (thermal  desorption).  In  flash  desorption  one  prepares  a  gas- 
solid  system  in  equilibrium,  and  the  barrier  is  surmounted  due  to  rapid 
heating  of  the  solid.  Isothermal  desorption  corresponds  to  suddenly  re¬ 
ducing  the  gas  pressure  while  keeping  the  solid  at  a  constant  temperature. 
Although  desorption  is  commonly  induced  by  vibrational  excitation,  it  can 
be  accomplished  by  electronic  excitation.  We  concentrate  here  on  the 
theory  of  the  former,  with  a  reference  to  electronic  excitation  at  the 
end  of  this  section. 

The  kinetics  of  desorption  is  conventionally  interpreted  in  terms 
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of  the  transition-state  theory.  Here  one  assumes  desorption  to 

proceed  via  an  energetically  activated  state  which  exists  at  the  top  of 
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a  potential  energy  barrier.  The  theory  yields  the  standard  Frenkel' 
Arrhenius  formula  for  the  rate  of  desorDtion, 


k  =  kQ  e'D/kBT, 


where  D  is  the  activation  energy,  kg  is  the  Boltzmann  constant,  T  is  the 

temperature  (the  final  temperature  of  the  solid  for  the  case  of  flash 

desorption),  and  the  preex^onenti al  factor  kQ  ranges  typically  from 
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10  /sec  to  10  /sec  depending  i  the  system.  Although  it  is  often 
assumed  as  a  first  approximation  that  kQ  and  D  are  constants,  they  in 
fact  vary  with  the  surface  coverage  due  to  interaction  b-tween  adatoms. 
Furthermore,  the  Frenkel -Arrhenius  formula  is  found  to  be  valid  only 
over  a  limited  range  of  temperature.  The  ultimate  goal  left  for  theorists 
to  achieve  to  find  the  rule  that  governs  the  dependence  of  the  desorp¬ 
tion  rate  on  temperature,  surface  coverage  and  other  possible  factors 
such  as  interaction  anharmoni city  and  lattice  dynamics.  This  has  to  be 
accompanied  by  an  accurate  determination  of  D  and  kQ  for  a  given  system. 
Dynamical  aspects  such  as  energy  and  angular  distributions  of  desorbed 
particles  can  also  be  studied.  In  the  following,  we  briefly  review  some 
selected  theoretical  progress  made  during  the  last  several  years. 

A.  Classical  Theory 

Perhaps  the  most  straightforward  way  of  computing  the  desorption 
rate  within  the  framework  of  classical  mechanics  is  to  employ  the  classical 
trajectory  method  widely  used  in  molecular  dynamics  calculations,  i.e., 
to  numerically  integrate  Hamilton's  equations  for  particles  in  the  system 
with  an  appropriate  distribution  of  initial  positions  and  momenta 
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characterized  by  a  given  temperature.  The  immediate  difficulty  of  this 
method  lies  in  the  fact  that  desorption  is  a  many-body  problem,  and 
consequently  the  number  of  coupled  equations  to  be  integrated  is  much 
too  large  for  the  method  to  be  of  any  use.  This  difficulty,  however, 
can  be  alleviated  by  the  GLE  approach  as  described  in  Sec.  II,  in  which 
only  a  few  surface  atoms  in  the  immediate  neighborhood  of  the  adatom  in 
question  are  given  a  full  consideration. 

An  early  attempt  at  a  kinetic  description  of  desorption  in  the 
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spirit  of  the  GLE  approach  was  made  by  Kramers.  Grimmelmann,  Tully  and 
45 

Helfand  applied  the  GLE  approach  to  compute  the  desorption  rate  of 
xenon  atom  from  the  (111)  surface  of  a  platinum  crystal.  The  primary  zone 
was  considered  to  contain  the  Xe  atom  and  four  nearby  Pt  surface  atoms 
whose  motion  was  followed  in  full  detail.  If  the  Xe  atom  moves  outside 
the  initially  designated  primary  zone,  original  primary  atoms  were 
switched  off  and  new  ones  switched  on.  Desorption  was  considered  to 
occur  when  the  separation  between  Xe  and  the  surface  is  greater  than 

o  o 

10A.  The  adatom-surface  potential  is  sufficiently  weak  beyond  10A 
that  particles  entering  this  region  may  be  assumed  to  be  desorbed.  The 
desorption  rate,  however,  was  reported  only  at  a  few  high  temperature 
values  (>500°K).  The  difficulty  is  that,  even  with  the  GLE  approach,  the 
computation  of  the  desorption  rate  requires  prohibitively  large  amounts 
computer  time.  This  is  because  desorption  is  fundamentally  an  infrequent 
event.  Although  trajectories  leading  to  desorption  must  pass  over  the 
top  of  the  potential  energy  barrier,  this  region  in  most  cases  is  re¬ 
presented  by  a  "bottleneck"  in  phase  space.  Trajectories  spend  most  of 
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the  time  wandering  in  the  low-energy  region  of  phase  space  and  almost 
never  reach  the  bottleneck  region.  Consequently,  desorption  is  a  slow, 
infrequent  event,  especially  at  low  temperatures.  It  therefore  is 
apparent  that  some  technique  of  efficiently  simulating  infrequent  events 
needs  to  be  employed  in  order  to  successfully  describe  desorption  processes 
within  the  framework  of  the  classical  Langevin  approach. 

One  such  technique,  a  compensating  potential  method,  was  proposed 
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by  Grinmelman,  Tully  and  Helfand.  In  this  method,  the  interaction 
potential  V,  which  in  general  depends  on  the  coordinates  of  the  adatom 
and  surface  atoms,  is  separated  into  two  terms,  V  =  W  +  U,  in  such  a  way 
that  U  depends  only  on  a  single  coordinate,  the  separation  between  the 
adatom  and  the  surface.  It  then  is  possible  to  rearrange  the  oroblem 
in  such  a  way  that  trajectories  subject  only  to  the  potential  W  need  to 
be  computed  for  the  evaluation  of  the  desorption  rate.  If  one  chooses  U 
so  that  it  contains  a  potential  minimum  of  approximately  the  same  depth 
as  V,  then  W  contains  only  a  shallow  minimum.  The  compensating  potential 
method  therefore  improves  the  efficiency  of  sampling  the  bottleneck  region 
by  effectively  replacing  the  deep  minimum  in  the  interaction  potential  by 
a  shallow  minimum.  With  this  technique,  the  desorption  rate  of  a  xenon 
atom  from  a  platinum  surface  was  computed  down  to  the  temperature  T  *  100°K. 
The  computed  rates  are  seen  to  be  consistent  with  the  Frenkel -Arrhenius 
formula  over  the  temperature  range  100°K  <  T  <  5C;’,‘..  However,  as  the 
temperature  is  increased  beyond  500°K,  the  rates  show  deviation  from  this 
simple  formula. 

Another  way  to  describe  infrequent  events  is  the  pseudo -dynamical 
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Monte  Carlo  sampling  technique.  Adams  and  Doll46,47  have  incorporated 
this  technique  into  the  classical  transition-state  theory  to  compute 
desorption  rates.  This  involves  finding  the  average  rate  at  which  parti¬ 
cles,  thermally  distributed  at  a  given  temperature,  enter  the  "desorption 
region"  defined  by  s  >_q,  where  s  is  the  adatom-surface  separation  and 
q  is  an  appropriately  chosen  separation  beyond  which  the  adatom-surface 
potential  is  sufficiently  weak.  Since  desorption  is  an  infrequent  event, 
it  is  not  practical  to  numerically  integrate  Hamilton's  equations  in 
order  to  generate  the  function  s  =  s(t).  The  Monte  Carlo  sampling 
technique,  however,  alleviates  this  difficulty  by  allowing  one  to  con¬ 
sider  "pseudo-dynamics"  of  adatoms.  The  technique  provides  an  efficient 
means  of  generating  averages  over  a  distribution  characterized  by  a 
given  temperature  and  is  most  suitable  for  the  description  of  desorption 
processes.  With  this  technique  it  is  possible  to  compute  desorption 

rates,  for  example,  for  a  wide  range  of  surface  coverages.  The  Monte 

46 

Carlo  transition-state  method  has  been  applied  to  desorption  of  a 
helium  atom  from  the  (111)  surface  of  xenon  with  the  He-Xe  interaction 
as  well  as  He-He  interaction  assumed  to  be  given  by  the  Lennard-Jones 
6-12  potential.  The  computed  rates  reveal  a  clear  dependence  on  the 
surface  coverage.  The  rate  was  found  to  be  an  increasing  function  of 
the  coverage.  For  example,  in  passing  from  the  infinite  dilute  limit 
to  a  monolayer  coverage,  the  rate  was  almost  doubled.  This  is  expected 
because  adatom-adatom  interactions  are  mainly  repulsive.  At  high  coverages, 
adatoms  tend  to  force  neighboring  atoms  off  the  surface,  leading  to  an 
increase  in  the  desorption  rate. 
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Another  interesting  feature  revealed  by  the  Monte  Carlo  calculation 
is  that  the  temperature  dependence  of  the  rate  shows  the  existence  of  a 
transition  temperature  Tc  =  30  ^  40°K  for  He-Xe(lll),  i.e.,  the  plot 
of  ink  vs.  T  consists  of  a  line  whose  slope  changes  at  T  =  T£.  This 
indicates  a  change  in  the  activation  energy  D  and  the  preexponential 
factor  at  the  transition  temperature. 

Only  recently  have  powerful  numerical  techniques  such  as  the  ones 
discussed  above  been  introduced  to  describe  desorption  processes,  thus 
making  possible  relatively  accurate  determination  of  desorption  rates  over 
a  wide  range  of  temperatures.  It  seems,  however,  that  many  of  the  dynamical 
details  of  desorption  still  remain  unanswered.  It  is  hoped  that  the 
development  of  a  comprehensive  theory  at  a  fundamental  level  as  well  as 
of  numerical  techniques  along  the  line  discussed  above  will  lead  to  a 
complete  understanding  of  a  desorption  process. 

B.  Quantum  Mechanical  Theory 

Quantum  mechanically,  thermal  desorption  is  pictured  as  a  transition 
of  an  adatom  from  a  bound  state  to  a  continuum  state,  with  the  required 
energy  supplied  by  the  phonon  bath  of  the  solid.  The  kinetics  of  desorp¬ 
tion  is  essentially  determined  by  a  series  of  transition  matrix  elements, 

Hmn  =  <mlH'ln>»  which  describes  a  transition  from  the  m-th  state  to  the 
n-th  state  (H'  is  the  interaction  Hamiltonian  for  the  adatom-surface). 

The  success  of  a  quantum  mechanical  theory  therefore  depends  largely  on 
the  capability  of  obtaining  accurate  eigenfunctions  to  evaluate  the  transi¬ 
tion  matrix  elements. 
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An  early  quantum  mechanical  theory  of  desorption  is  due  to  Lennard- 
Jones,  Strachan  and  Devonshire.4®*49  They  chose  a  Morse  potential  to 
represent  an  adatom-surface  atom  interaction  and  evaluated  the  transition 
matrix  elements  accordingly.  The  Morse  potential  seems  particularly 
useful  here  because  the  corresponding  eigenfunctions  are  known  and  the 
matrix  elements  can  be  analytically  evaluated.  In  applying  their  formalism 
to  desorption,  however,  they  have  limited  their  consideration  to  the  case 
where  the  activation  energy  is  sufficiently  low  that  a  direct  transition 
from  the  initial  bound  level  to  a  continuum  level  can  be  mediated  by  a 
single  phonon.  The  theory  therefore  is  applicable  to  only  a  limited  num¬ 
ber  of  systems,  with  a  very  low  activation  energy.  In  addition,  their 
theory  is  essentially  one  dimensional  in  that  the  motion  of  the  adatom 
or  surface  atoms  parallel  to  the  surface  olays  no  role. 

Bendow  and  Ying50  have  attempted  a  generalization  of  this  early 
theory  by  considering  multiphonon  desorption  in  three  dimensions.  From 
a  practical  viewpoint,  the  generalization  implies  a  greater  number  of 
transition  matrix  elements  to  be  evaluated.  Furthermore,  each  element 
is  now  more  complicated  because  the  eigenfunctions  depend  not  only  on 
the  adatom  surface  separation  but  also  on  the  coordinates  parallel  to 
the  surface  which  reflect  the  periodicity  of  the  lattice  structure. 

Applying  their  theory  to  neon  desorbed  from  xenon-covered  graphite, 
they  found  that  the  preexponential  factor  kQ  is  unusually  small, 

5 

k0  10  /sec.  It  was  suggested  that  this  may  be  due  to  the  inherent 
difference  between  one-dimensional  and  three-dimensional  models. 
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Goodman  and  Romero  reported  calculations  of  desorption  rates  for 
the  systems  H  on  C  and  He  on  Ar,  Kr,  Xe  assuming  that  desorption  occurs 
via  a  single-phonon  transition.  They  used  a  three-dimensional  quantum 
mechanical  approach  and  assumed  that  the  adatom-surface  interaction  is 
given  by  a  Morse  potential  with  appropriately  chosen  parameters  for  each 
system.  In  all  cases  the  rate  is  seen  to  follow  the  Frenkel -Arrhenius 
formula  with  kQ  ranging  from  10^/sec  to  10^/sec.  Calculations  have 
also  been  carried  out  with  one-dimensional  models.  In  contrast  to  the 
suggestion  by  Bendow  and  Ying,  however,  the  one-dimensional  and  three- 
dimensional  models  yielded  essentially  the  same  desorDtion  rates. 

A  simple  model  of  multiphonon  desorption  was  put  forth  by  Garrison, 
52 

Diestler  and  Adelman  in  which  a  one-dimensional  truncated  harmonic 
potential  is  used  to  describe  the  atom-surface  interaction,  and  eigen¬ 
functions  are  assumed  still  to  be  given  by  those  of  the  untruncated 
harmonic  potential.  The  transition  matrix  elements  can  then  be  easily 
evaluated.  The  rate  constant  was  simply  taken  as  the  reciprocal  of  the 
time  required  to  reach  the  continuum.  This  model  was  applied  to  desorpti 
of  a  xenon  atom  from  a  tungsten  surface.  The  ca’culated  desorption  rate 
exhibits  Frenkel -Arrhenius  behavior  and  depends  on  lattice  dynamical 
properties;  in  particular,  it  shows  an  inverse  cube  dependence  on  the 
Debye  temperature.  Although  the  desorption  rate  is  expected  to  depend 
upon  lattice  dynamics,  unrealistic  assumptions  made  in  this  model 
(truncated  harmonic  potential,  harmonic  wavefunctions)  limit  the  reliabi¬ 
lity  of  calculated  data. 
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Gortel,  Kreuzer  and  Teshima  have  developed  a  quantum  mechanical 
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theory  based  on  the  rate  equation 


dn^t) 
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where  n..  is  the  time-dependent  occupation  of  the  i-th  bound  state,  R^ 
is  the  rate  of  transition  from  the  i-th  to  j-th  bound  state,  and  Rc- 
is  the  transition  rate  from  the  i-th  bound  state  into  the  continuum. 

The  theory  applies  mainly  to  systems  with  low  coverage,  because  non¬ 
linear  terms  that  account  for  the  adatom-adatom  interaction  are  not 
included  in  Eq.  (88).  It  can  be  shown  that  the  desorption  rate  is  the 
smallest  eigenvalue  of  the  matrix  R^.  The  problem  then  is  reduced  to 
finding  the  R. .,  which  in  turn  are  determined  by  the  matrix  elements 

*  J 

H^.  The  matrix  elements  were  evaluated  exactly,  assuming  a  Morse 

potential  between  the  adatom  and  surface.  Calculations  showed  that 

desorption  in  weakly  coupled  systems  with  many  bound  states  proceeds 

predominantly  through  one-phonon  cascades.  Isothermal  desorption  rates 

were  calculated  for  He-LiF,  He-NaF,  He-graphite,  H-NaCl,  He-Ar  and  Xe-W. 

All  the  calculated  rates  were  found  to  behave  according  to  the  Frenkel  - 

Arrhenius  formula  with  activation  energy  approximately  equal  to  the  energy 

54 

of  the  lowest  bound  state.  Gortel,  Kreuzer,  Teshima  and  Turski  showed 
that  a  further  simplification  results  if  the  adatom-surface  potential 
allows  a  large  number  of  bound  states  such  as  in  the  Xe-W  system.  In 
this  case,  the  system  may  be  viewed  as  having  a  quasi -continuum  instead 
of  many  discrete  bound  states.  The  system  may  then  be  considered  to 
perform  a  random  walk  through  this  quasi -continuum  rather  than  cascading 


through  bound  states.  The  rate  equation,  Eq.  (88),  can  then  be  cast  into 
a  form  in  the  continuum  limit,  and  from  this  the  Fokker-Planck  equation 
derived.  This  equation  was  shown  to  lead  to  a  relatively  simple  expression 
for  the  desorption  rate  for  the  case  of  a  weakly-coupled  system. 

In  computing  the  desorption  rate,  one  needs  to  carry  out  thermal 

averaging  of  the  square  of  the  transition  matrix  elements  over  Dhonon 

modes.  At  this  stage,  virtually  all  the  theories  that  describe  desorption 

have  used  a  bulk  Debye  model  for  the  solid,  i.e.,  it  was  assumed  that 

phonons  that  induce  desorption  of  adatoms  from  a  surface  are  simply  those 

of  an  infinite  solid.  However,  surface  phonon  modes  may  significantly 

differ  from  those  of  solid,  and  therefore  a  proper  account  of  such  modes 

should  be  given.  This  problem  was  undertaken  by  Goldys,  Gortel and 
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Kreuzer,  who  found  that  surface  phonons  contribute  to  desorption  two 
to  three  times  as  much  as  bulk  phonons.  From  this,  an  effective  surface 
Debye  temperature  can  be  deduced,  which  turns  out  to  be  about  0  65  to 
0.8  of  the  bulk  Debye  temperature. 

Despite  progress  made  during  recent  years,  many  details  of  desorption 
are  still  uncharacterized.  One  of  the  important  aspects  of  desorption  not 
yet  totally  answered  is  the  dependence  of  the  rate  upon  surface  coverage. 
Adatom-adatom  interactions  are  often  neglected,  which  makes  it  difficult 
to  assess  the  effect  of  surface  coverage  on  desorption.  The  relation 
between  the  desorption  rate  and  the  adatom-surface  interaction  potential,  e.g. 
the  effect  of  interaction  anharmonicity  on  the  desorption  rate,  still  seems 
unclear.  Although  calculations  of  desorption  rates  for  different  Interaction 
potentials  exist,  it  appears  that  no  systematic  treatment  of  this  problem 
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has  been  given.  Nevertheless,  some  accurate  calculations  of  the  rates 
at  low  coverages  over  a  wide  range  of  temperatures  have  begun  to  appear 
in  recent  years. 

Finally,  we  mention  that  desorption  can  also  be  induced  by  electronic 
transitions,  with  the  use  of,  say,  an  electron  or  photon  beam.  In  order 
for  desorption  to  be  induced  by  electronic  transitions,  an  electron  must 
be  excited  to  a  state  with  a  repulsive  potential  energy  surface.  The 
molecules  or  ions  will  then  be  repelled  and  rejected  from  the  surface. 

The  main  issue  here  is  not  the  motion  of  electrons  but  the  nuclear  dynamics 
that  follows  the  electronic  excitation.  For  a  detailed  treatment  of  de- 
sorption  induced  by  electronic  transitions,  we  refer  to  Gadzuk.  We 
shall  also  consider  this  problem  in  Sec.  VI  in  terms  of  surface  state 
excitation  within  the  substrate. 
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VI.  LASER-STIMULATED  RATE  PROCESSES 

57  5ft 

Both  the  phonon  and  electron  spectra  of  a  solid  with  a  surface 
have  energy  levels  distinct  from  the  bulk  levels  of  an  infinite  solid. 
These  levels  correspond  to  vibrations  or  electronic  charge  density 
localized  in  the  vicinity  of  the  surface.  The  introduction  of  an 
adspecies  on  this  surface  will  also  produce  phonon  and  electron  energy 
levels  related  to  the  adsorptive  bond. 

Using  a  laser,  the  phonon  mode  of  the  adsorptive  bond  could  be 
excited  in  a  selective  manner  in  order  to  induce  desorption,  migration 
and  bond-breaking  within  the  adspecies.  A  thermal  non-selecti ve  excita¬ 
tion  of  the  surface  phonon  modes  by  the  laser  could  also  be  used  to 
stimulate  these  dynamical  processes,  in  addition  to  serving  as  a  means 
of  annealing.  In  Part  A  of  this  section,  we  shall  discuss  the  laser 
excitation  of  these  adsoecies-surface  phonon  modes  in  detail.  After  this 
discussion,  we  shall  examine  in  Part  B  the  effects  of  laser  radiation  on 
the  excitation  of  surface  electrons  on  a  semiconductor.  It  will  be  shown 
that  such  excitation  can  have  a  significant  influence  on  the  adsoecies- 
surface  interaction  potential.  Part  C  will  demonstrate  the  importance 
of  including  both  electrons  and  phonons  in  the  laser  excitation  of  a 
metal  surface.  Finally,  in  Part  D  we  will  discuss  the  effects  on  charge 
transfer  between  the  surface  and  an  adspecies  produced  by  the  altered 
surface  charge  caused  by  laser  excitation. 

A.  Adspecies-Surface  Vibrational  Modes 

The  vibrational  spectrum  of  an  adspecies-surface  system  can  be 


divided  into  a  number  of  distinct  modes.  The  active  modes  of  the 
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adsorptive  bond  along  with  intramolecular  vibrations  of  the  adspecies 
are  different  from  the  bath  modes  of  the  solid.  The  bath  modes  in  turn 
can  be  viewed  as  having  components  due  to  surface  and  bulk  vibrations 
of  the  solid.  Fig.  1  is  a  schematic  representation  of  such  a  spectrum. 

An  active  mode  (A)  is  indicated  by  the  highest  energy  peak;  the  high- 
energy  bath  modes  (B)  would  be  those  due  to  surface  vibrations;  and  the 
low-energy  modes  (C)  are  due  to  the  internal  vibrations  of  the  solid. 

Laser  irradiation  of  the  solid  surface  with  a  frequency  tuned  to 
the  active  mode  will  lead  to  energy  transfer  to  the  adsorptive  bond. 

This  can  produce  selective  bond  breaking  between  the  adspecies  and  sur¬ 
face  or  within  the  adspecies.  However,  to  accomplish  this  bond  breaking, 
the  laser  must  pump  energy  into  the  active  mode  faster  than  the  vibration 
is  damped  by  the  coupling  between  the  active  and  bath  modes.  This 
vibrational  relaxation  is  characterized  by  both  an  energy  (T^)  and  phase 
(T^)  process.  In  addition  to  this  anharmonic  coupling  to  the  bath  phonon 
modes,  the  T1  and  Tj  processes  are  also  effected  by  vibration-induced 
migration  of  the  adspecies,  substrate-induced  thermal  fluctuations  of 
the  effective  dipole  of  the  adspecies,  and  charge  transfer  between  the 
surface  and  the  adspecies. 
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Both  a  classical  model  based  on  the  GLE  approach  and  a  quantum 
model  based60  on  the  Heisenberg-Markov  approximation  have  been  used  to 
describe  the  relaxation  process  of  these  adspecies-surface  dynamics. 

The  results  have  shown  that  for  a  given  anharmonicity  of  the  effective 
adsorptive  potential,  given  detuning  between  the  laser  frequency  and  the 
active  mode  frequency,  and  given  dephasing  lifetime,  then  selective 
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excitation  is  confined  to  a  small  range  of  values  of  the  danping  factor 

y  1  which  is  associated  with  T-j  processes.  These  selective  dynamics,  on 

the  other  hand,  are  rather  insensitive  to  T£  processes. 

Fig.  2  shows  the  effects  of  the  energy  relaxation  rate  Yj  on  the 

amount  of  selective  excitation  of  the  active  mode  at  the  expense  of  the 
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bath  modes.  As  can  be  seen  by  the  lower  curve  (E),  there  is  a  threshold 

for  Yj  beyond  which  essentially  all  selectivity  disappears.  Also,  low 

selectivity  for  short  irradiation  times  followed  by  enhanced  selectivity 

for  large  exposure  times  points  to  important  feedback  mechanisms  as  well 

as  an  optimum  irradiation  time  for  maximum  selectivity. 

Despite  its  importance,  values  of  yj  are  not  easily  obtainable  from 

either  theory  or  experiment.  A  mostly  first-principles  formulation  for 
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y ^ ,  however,  has  recently  been  developed.  The  relaxation  rate  is  given 
by  the  Fermi  golden  rule  for  a  perturbation  which  is  the  difference  in 
the  potential  due  to  the  lattice  atoms  at  their  instantaneous  and  equili¬ 
brium  positions.  This  rate  reduces  to  a  sum  of  products  of  matrix  elements 
involving  only  adatom  properties  and  time  Fourier  transforms  of  correlation 
functions  for  different  lattice  displacements  at  different  times. 

In  addition  to  the  selective  excitation  of  phonon  modes,  laser- 
stimulated  surface  processes  of  a  non-selecti ve,  thermal  nature  can  be 
important  in  such  areas  as  annealing.  The  techniques  for  describing 

selective  excitations  can  be  combined  with  a  diffusion  equation  to  treat 

6  ^ 

overall  energy  transfer  processes.  Both  selective  and  non-selecti ve 
energy  transfer  via  laser  excitation  of  the  adspecies-surface  phonon 
modes  is  seen  to  be  important.  The  next  step  in  the  theory  of  laser 
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induced  surface  dynamics  would  involve  the  inclusion  of  electronic 

degrees  of  freedom. 

B.  Charging  a  Semiconductor  Surface 

Laser  excitation  of  the  electronic  degrees  of  freedom  in  the  solid 

can  stimulate  rate  processes  by  altering  the  charge  in  the  vicinity  of 

a  surface.  Such  alteration  of  charge  can  induce  desorption  of  an 

adspecies.  Experiments  with  synchrotron  radiation  on  metal  surfaces  has 

64 

demonstrated  this  phenomena  for  a  variety  of  adspecies.  It  is  reason¬ 
able  to  assume  that  reversing  the  charge  on  the  surface  could  also 
enhance  adsorption  of  the  same  species. 

To  analy2e  this  process,  a  truncated  one-dimensional  chain  has 
65 

been  employed.  Within  this  model  and  the  nearly-free-electron  approxi¬ 
mation,  the  electronic  energies  are  given  by 

Ek  =  l  {[kZ  +  (k  -  g)2]  ±  /  [k2  -  (k  -  g)2]2  +  4E2  }  , 

where  k  is  the  wave  vector  of  the  electron,  g  is  the  reciprocal  lattice 
constant,  and  Eg  is  the  band  gap.  The  plus  branch  represents  the  con¬ 
duction  energy  band  and  the  minus  branch  corresponds  to  the  valence 
band.  Surface  states  are  obtained  by  the  use  of  analytical  continua¬ 
tion,  i.e.,  the  energy  can  also  have  real  values  for  k  =  g/2  +  i*. 

Since  1/2*  is  the  effective  charge  depth  of  these  states,  they  are  most 
like  surface  states  for  large  values  of  <.  The  band  structure  for 
this  system  is  given  in  Fig.  3.  The  band  extending  out  of  the  figure 
represents  the  possible  surface  states. 

If  a  laser  is  now  shown  on  the  linear  chain,  we  can  transfer 
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electrons  from  the  valence  band  to  the  surface  states.  Considering  the 

Bloch  nature  of  our  bulk  states  and  the  damped  Bloch  nature  of  our  surface 

states,  a  simple  analysis  of  the  transition  matrix  will  lead  to  the 

conclusion  that  the  transfer  is  favored  for  bulk  states  at  the  band  edge. 

Using  this,  we  have  calculated  the  absorption  cross  section  for  laser 

excitation  of  a  silicon  surface,  and  our  results  are  depicted  in  Fig.  4. 

The  glitch  near  the  center  corresponds  to  the  maximum  value  of  <  and  is 

due  to  the  branch  point  in  the  energy.  Since  k  is  a  maximum  near  this 

region,  to  achieve  the  largest  surface  charge  density,  surface  states 

near  this  center  should  be  excited.  We  can  see  from  the  plot  that  quite 

°2 

large  cross  sections  (2-3A  )  can  be  obtained  in  this  region  at  low  radia- 

2 

tion  densities,  (<10  W/cm  ).  Thus  a  low-power  laser  is  a  very  effective 
controller  of  surface  charge. 

To  determine  the  effect  of  this  charge  on  an  adspecies,  we  have 
looked  at  the  interaction  of  a  charged  adspecies  with  the  surface.66 
The  exact  form  of  the  screening  function^7  at  the  surface  can  be  extremely 
complicated;  therefore,  as  a  first  approximation,  we  will  assume  that  the 
screening  is  represented  by  an  exponential  damping  (Thomas -Fermi  screening  ). 
Our  results  for  the  change  in  the  interaction  for  a  number  of  different 
surface  state  excitations  are  plotted  in  Fig.  5.  It  should  be  noted  that 
states  near  the  center  of  the  band  have  a  great  deal  more  interaction 
as  would  be  expected  because  of  the  smaller  charge  depth.  As  can  be  seen 
from  the  figure,  both  the  magnitude  and  the  range  of  the  interaction  are 
increased  by  the  laser-charged  surface.  The  attraction  of  the  ion  for 
the  surface  will  thus  be  greatly  enhanced  if  the  ion  is  positive  and 
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diminished  if  the  ion  is  negative.  A  similar  analysis  could  be  done 
for  polar  species,  with  one  pole  being  attracted  and  the  other  repelled 
Furthermore,  a  consideration  of  a  three-dimensional  semiconductor  would 
also  introduce  occupied  surface  states.  Thus  raising  the  possibility 
of  exciting  holes  into  the  surface,  subsequently  making  it  more  positive, 
and  producing  the  reverse  effect  of  electronic  excitation. 

C.  Electron-Phonon  Coupling  on  Metals 

In  the  preceding  Part  B,  we  looked  at  the  role  of  increasing  sur¬ 
face  charge  in  semiconductors  and  its  effect  on  surface  dynamics.  Since 
metals  also  have  surface  states,  it  would  be  expected  that  similar  charge 
transfer  induced  with  a  laser  would  be  possible  in  metallic  systems. 

Again  using  a  truncated  one  dimensional  modeled  the  nearly-free- 
electron  approximation, 69  the  band  structure  for  the  metal  will  be  found 
qualitatively  similar  to  that  of  a  semiconductor  (see  Fig.  3).  However, 
whereas  the  lower  band  is  populated  to  the  top  in  a  semiconductor,  this 
band  in  a  metal  would  be  only  partially  filled.  For  example,  for  the 
metal  sodium  the  top  of  the  band  is  0.7  eV  above  the  Fermi  energy. 
Consequently,  to  excite  surface  states  in  a  metal,  we  must  first  overcome 
this  energy  deficit  by  application  of  a  higher  frequency  laser  than  needed 
for  semiconductors. 

In  a  calculation  of  the  cross  section,  we  found  that  our  selection 
rule  requiring  conservation  of  the  real  part  of  the  crystal  momentum 
still  holds  for  metallic  system.  Since  the  top  of  the  band  is  unpopulated, 
we  would  therefore  expect  zero  transitions  in  the  metal.  However,  the 
introduction  of  phonons  can  overcome  this  problem.  It  should  be  noted 
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that  neither  laser  photons  or  lattice  phonons  are  alone  able  to  populate 
the  surface  states.  The  former  will  not  supply  sufficient  crystal 
momentum,  and  the  latter  will  not  supply  sufficient  energy.  Consequently, 
the  absorption  cross  section  will  be  second  order  in  nature  with  a 
mixing  of  phonon  and  photon  excitations.  The  dominant  pathway  would  be 
a  phonon  first  scattering  the  electron  from  the  top  of  the  Fermi  energy 
to  the  top  of  the  band  and  then  a  laser  photon  exciting  this  electron 
to  the  desired  surface  state. 

The  mathematics  involved  in  the  calculation  of  the  cross  section 
will  be  similar  to  that  of  the  semiconductor.  Because  of  the  phonons, 
however,  an  averaging  over  the  populated  phonon  states  must  be  performed. 
This  will  give  a  factor  depending  on  the  phonon  population  multi Dlying 
our  first  order  results: 

The  first-order  cross  section  will  be  comparable  to  that  of  silicon, 

given  in  Fig.  4.  For  sodium,  the  scaling  factor  S  is  about  10~4. 

Consequently,  the  statements  that  applied  for  semiconductors  would  also 

be  applicable  to  metals  if  the  power  density  of  our  laser  were  104 

times  as  large.  Since  large  cross  sections  were  observed  in  the  semi- 

2 

conductor,  only  a  low-power  laser  (1-10  W/cm  )  was  necessary.  Therefore, 

2 

only  a  moderate  power  laser  (10-100  kW/cm  )  would  be  needed  to  induce 
surface  charge  transfer  in  a  metal.  Consequently,  processes  such  as 
enhanced  or  diminished  adsorption  or  desorption  of  charged  or  polar 
adspecies  would  also  occur  with  similar  results  for  metal  or  semiconductor 
surfaces. 
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D.  Gas -Surface  Charge  Transfer 

If  a  laser  is  used  to  charge  a  surface,  we  would  expect  the  transfer 
of  charge  between  the  solid  and  an  adspecies  or  gas  molecule  above  the 
surface  to  be  greatly  enhanced.^  Since  charge  transfer  is  the  essence 
of  chemistry,  surface  catalysis  would  be  greatly  effected  by  this 
phenomenon. 

If  we  consider  a  positive  ion  impinging  on  a  semiconductor  surface, 
the  wave  function  of  the  entire  system  could  be  approximated  as  a  product 
of  the  wave  function  of  the  crystal  and  the  ion.  On  leaving  the  surface, 
the  ion  could  be  neutralized  by  the  transfer  of  charge  from  the  surface. 
The  final  state  would  thus  be  a  product  of  the  wave  functions  of  the 
neutral  gas  molecule  and  the  charged  substrate.  If  the  surface  of  the 
semiconductor  is  first  charged  with  a  laser,  the  final  state  with  the 
neutral  gas  molecule  would  be  very  favored.  Fig.  6  illustrates  the  range 
over  which  charge  transfer  with  the  surface  states  can  occur.  The  shaded 
curve  in  the  figure  is  due  to  the  thickness  of  the  surface  band  as  shown 
in  Fig.  3.  Thus,  if  the  ion  is  on  its  initial  potential  curve,  Hjj,  the 
charge  neutralization  process  will  be  enhanced  where  this  curve  intersects 
the  neutral  gas-charged  crystal  interaction  illustrated  by  the  shaded 
curves.  A  mathematical  analysis  using  first-order  time-dependent  pertur¬ 
bation  theory  will  show  that  the  neutralization  rate  of  the  ion  is  larger 
by  a  factor  of  j 2<L j  for  a  surface  electron  in  state  <  than  for  a  bulk 
electron  (L  is  the  length  of  the  solid). 
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Figure  Captions 

Fig.  1.  Schematic  diagram  of  the  phonon  density  of  states.  A  is  the 
active  mode,  and  B  and  C  are  bath  modes. 

Fig.  2.  Average  selective  excitation  number  <N^>  as  a  function  of  time  t 
for  various  values  of  the  energy  relaxation  rate  y-j  for  a  model 
adspecies-surface  system:  (A)  10^s”^;  (B)  103s_1;  (C)  2  *  lO^s”1; 

(D)  6  x  loV1;  and  (E)  loV1. 

Fig.  3.  Dispersion  relationship  in  complex  crystal  momentum  space  (k  +  i<) 
for  a  finite  linear  chain.  The  valence,  surface  and  conduction 
bands  are  labeled  V,  S  and  C,  respecti vely . 

°2 

Fig.  4.  Absorption  cross  section  for  surface  states,  o,  in  A  versus 

the  frequency  of  the  exciting  laser  radiation. 

Fig.  5.  Magnitude  of  the  surface  interaction  potential  (in  mi  1 1 ihartrees )  at 

various  distances  from  the  surface.  The  solid  line  represents  the 

system  with  excited  state  <  £  Eg/g ;  the  dashed  line,  k  =  -  0.5  Eq/a; 

and  the  dotted  line,  k  =  -  0.1  Eg/g,  all  in  the  lower  energy  branch. 

Fig.  6.  Schematic  diagram  of  potential  energy  surfaces  Hjj  and  Hpp  for 

initial  and  final  configurations.  The  transfer  of  a  surface 

electron  to  the  incident  ion  occurs  at  time  t-j  <_  t  <_  t2«  For 

neutralization  of  He  at  Si(lll)  into  the  2  S  state  of  He,  A  n-  0.07  eV 

and  Eg  ^  1.1  eV.  Shown  in  the  inset  is  the  square  of  the  normaliza- 
2 

tion  constant,  N^.  It  assumes  a  maximum  value  when  the  energy  E  is 
at  the  center  of  the  gap  and  approaches  zero  at  both  ends  of  the 
gap.  This  indicates  that  surface  electrons  with  energy  near  the 
center  of  the  gap  are  most  highly  localized  near  the  surface. 
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